Transport properties of ID disordered models: a novel approach 
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A new method is developed for the study of transport properties of ID models with random 
potentials. It is based on an exact transformation that reduces discrete Schrodinger equation in the 
tight-binding model to a two-dimensional Hamiltonian map. This map describes the behavior of a 
classical linear oscillator under random parametric delta-kicks. We are interested in the statistical 
properties of the transmission coefficient Tl of a disordered sample of length L. In the ballistic regime 
we derive expressions for the mean value of the transmission coefficient Tl, its second moment and 
variance, that are more accurate than the existing ones. In the localized regime we analyze the global 
characteristics of In Ti,, and demonstrate that its distribution function approaches the Gaussian form 
if L — > oo. For any finite L there are deviations from the Gaussian law that originate from the subtle 
correlation effects between different trajectories of the Hamiltonian map. 



I. INTRODUCTION 



In recent years the interest to ID models with random 
potentials has been significantly increased. There are two 
reasons for this. First, it is expected that ID systems may 
elucidate the origin of the famous single parameter scal- 
ing (SPS) [1] for transport characteristics of disordered 
conductors. As was claimed in Refs. [2], a random char- 
acter of the fluctuations of the Lyapunov exponent for 
finite-length samples, that was originally used [3] to jus- 
tify SPS, is not correct. Specifically, it was shown that 
in the vicinity of the band edges the SPS does not hold 
[2]. This result is important not only from the theoreti- 
cal viewpoint, but also for the experiment (see discussion 
and references in [2]). 

Another reason of the growing interest to ID random 
models is due to recent results on the correlated dis- 
order. In early studies of the ID tight-binding model 
with a specific site potential (the so-called random dimcr 
model [4]), it was demonstrated a highly non-trivial role 
of short-range correlations. It was found that for discrete 
values of energy E that are determined by the model pa- 
rameters, a random dimer turns out to be fully transpar- 
ent (for other examples, see [5]). Practically, this leads 
to an emergence of a finite range of energies close to the 
resonant one, where the localization length of eigen- 
states is larger than the size L of finite samples. Since the 
number of such states is of the order of \fh, this effect 
was claimed to have practical importance in application 
to polymer chains. 

Although the localization length l^ diverges at dis- 
crete values of E, this model does not exhibit a mobility 
edge. However, the situation was found to be very dif- 
ferent for the case of long-range correlations [6]. As was 
shown in Ref. [7], one can construct such random (cor- 
related) potentials that result (for a weak disorder) in a 
band of a complete transparency. The position and the 



width of the window of transparency can be controlled 
by the form of the binary correlator of the weak random 
potential [8]. 

The role of long-range correlations has been studied 
in details for the tight- binding Anderson-type model [7], 
and for the Kronig-Penney model with randomly dis- 
tributed amplitudes [9] and positions of delta-peaks [10]. 
The results have been also extended to a single-mode 
waveguide with random surface profiles [11]. The pre- 
diction of the theory [7] has been verified experimen- 
tally [12], when studying the transport properties of a 
single-mode electromagnetic waveguide with point-like 
scatterers. The latter have been intentionally inserted 
into the waveguide in a way to provide a random poten- 
tial with slowly decaying binary correlator. Very recently 
[13] the existence of the mobility edges was predicted for 
waveguides with a finite number of propagating channels 
(quasi- ID system) with long-range correlations in ran- 
dom surface scattering potential. 

An appropriate tool to study the correlated disorder in 
ID Anderson and Kronig-Penney models is the so-called 
Hamiltonian map (HM) approach. The key point of 
this approach is a transformation that reduces a discrete 
ID Schrodinger equation to the classical two-dimensional 
Hamiltonian map. The properties of trajectories of this 
map are related to transport properties of a quantum 
model. The geometrical aspects of the HM approach turn 
out to be helpful in qualitative analysis as well as in de- 
riving analytical formulas. 

Originally, the HM approach was proposed in Ref. [14] 
in order to study extended states in the dimer model. 
It was demonstrated that with the help of this approach 
it is easy to understand the mechanism for the emer- 
gence of the extended states. Specifically, the conditions 
for resonant energies can be easily obtained not only for 
dimers, but also for N— mers (when blocks of the length 
N with the same site energy e appear randomly in a site 
potential). More important is that the expression for the 
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localization length loo in the vicinity of the band edges 
has been derived in a more general form, as compared to 
that found in earlier studies. 

The generalization of the HM approach to correlated 
ID-potentials [7,9,10] has led to an understanding of the 
role of long-range correlations for the emergence of ex- 
tended quantum states in random potentials. As a result, 
a new principle was proposed for filtering of stochastic 
and digital signals using random systems with correlated 
disorder [8]. 

So far, the HM approach was used for calculations of 
the localization length in infinite samples. A more serious 
problem arises when studying the transport properties of 
finite samples of length L. It is assumed [3] that sta- 
tistical properties of the transmission coefficient Tl are 
entirely determined by the finite-length Lyapunov expo- 
nent (FLLE) which fluctuates with different realizations 
of the random potential. According to this point, many 
studies of solid state models with random potentials are 
directly related to the analysis of statistical properties of 
the FLLE. 

In this paper we propose a new method based on the 
HM approach, that allows to obtain main transport char- 
acteristics in ID random models. To illustrate this ap- 
proach, we take the standard Anderson model with weak 
white-noise potential, and consider two limit cases of the 
ballistic and localized transport. We show that in this 
way one can relatively easily derive some of known re- 
sults, as well as obtain the new ones. 

II. THE HAMILTONIAN MAP APPROACH 

It is well known that the discrete ID Anderson model 
can be written in the following form of Schrodingcr equa- 
tion for stationary eigenstates ip n , 

tpn+i + ipn-i = {E + e n )ip n . (1) 

Here E is the energy of a specific eigenstate, and e n is the 
potential energy at site n. It is convenient to represent 
Eq.(l) in the form of two-dimensional Hamiltonian map 
[14,15], 

f x n +i \ _ ( cos it + y4„ sin it sin/i \ / x n \ , , 
\Pn+l ) V An cos A* _ sin M COSfl)\p n )' 

In this map the canonical variables x n = tp n and p n — 
(tpnCOSfi — VVi-i) / sin [i correspond to the position and 
momentum of a linear oscillator subjected to linear time- 
periodic delta- kicks. The amplitudes A n of the kicks 
are proportional to the site potential in Eq.(l), A n = 
—e n / sin/i, and the parameter it is determined by the 
energy of an eigenstate, 2 cos /i = E. 

One can see that the representation (2) is the Hamil- 
tonian version of the standard transfer matrix method. 
Indeed, starting from initial values Vo an d ip-i, one can 



compute ip n and ip n —i according to Eq.(l) or Eq.(2). The 
Lyapunov exponent A and, therefore, the localization 
length loo = A -1 is obtained in the limit n — > oo (sec 
below) . It turns out that the Hamiltonian representation 
(2) is more convenient than the standard one based on 
Eq.(l). This is due a possibility to introduce the clas- 
sical phase space in order to study the properties of a 
trajectory {x n ,p n }. 

Mathematically, the Anderson localization corre- 
sponds to the parametric instability of a linear oscilla- 
tor associated with the map (2), see details in [16] and 
some applications in [17]. One should stress that the 
unbounded classical trajectories of the map (2) do not 
correspond to the eigenstates of Eq.(l), however, they 
give a correct value of the localizaion length through the 
Lyapunov exponent. On the contrary, if an eigenstate 
ip n of Eq.(l) is extended (delocalizcd) in the infinite con- 
figuration space {n}, the corresponding trajectory of the 
map (2) is bounded in the phase space {x n ,p n }. In this 
case the trajectory specified by the value it has a direct 
correspondence to the eigenstate with energy E. There- 
fore, the structure of such eigenstates can be studied by 
analyzing the properties of the trajectories in the phase 
space. 

For our analysis it is convenient to represent the map 
of Eq. (2) in the action-angle variables (r, 9) . Using the 
standard transformation x = r sin# and p — rcos9, the 
map can be rewritten as follows, 

sin 9 n+1 = D^ 1 [sin (9 n - it) - A n sin 9 n sin /i] , (3) 
cos 9 n+ i = D^ 1 [cos {6 n - fi) + A n sin 9 n cos /j] 
where 



D n = \J 1 + A n sin 29 n + A\ sin 2 9 n . (4) 

Eqs.(3) and (4) allow one to represent the Lyapunov 
exponent A in terms of A n and 9 n (for energies not close 
to the band edges, see details in [15]), 

A = I- 1 = lim f V In ( rs±l) 

°° i-oo ti V r " ' (5) 

= i(ln(l + A„sin20„ + A^sin 2 9 n )) n . 

Here the angle brackets (...)„ stand for the average along 
a trajectory (over n). 

The relation (5) is valid for arbitrary potential e„, no 
matter, weak or strong, random or deterministic, pro- 
vided that the parameter \i corresponds to the value of 
E taken inside the energy spectrum. Note that for a weak 
potential the spectrum remains unperturbed (| E |< 2) 
in the lowest (Born) approximation. 
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From Eqs.(3) and (4) it follows that actually the small 
parameter is A n = —e n / sin rather than | e n |<C 1. 

Therefore, close to the edges of the energy band where 
| sin^| f« 0, the standard perturbation theory fails. That 
is why for energies close to the band edges one needs to 
use specific methods for calculation of the localization 
length (see [15] and references therein). 

If the energy is not close to the band edges, | E |= 2, 
or to the band center, E = 0, the standard perturba- 
tion theory with respect to A n is applicable. For a weak 
uncorrelated potential ("white noise") the distribution of 
9 n is homogeneous, V{9) = 1/2tt, i.e., the phases 9 n are 
independent of the potential e n . Therefore, instead of 
the average along the trajectory, one can perform an en- 
semble average over e n and 9 n independently The re- 
sult for loo can be obtained easily, by keeping the linear 
and quadratic terms in the expansion of the logarithm in 
Eq.(5), 



r 1 - 



! sin 2 n 



(6) 



Here and below the brackets (...) stand for the average 
over disorder. 

For the first time the result (6) was obtained by Thou- 
less [18] by another method. It is interesting to note that 
at the center of the energy band, E = 0, the correct 
expression for the localization length is slightly different 
from Eq.(6). At this point the Born approximation is 
invalid since the kinetic energy vanishes. 

Using the HM approach it is easy to see that at E = 
(where /i = 7r/2), classical trajectories reveal a mixture of 
the periodic rotation with period 4 (in number of kicks) 
around the origin p = x = 0, and a very slow diffusion in 
9 and r. As a result, the distribution function V(9) turns 
out to be slightly modulated over 9 by a periodic func- 
tion with period it/ 2 (see details and discussion in [15]). 
This leads to an anomalous contribution of the fourth 
Fourier harmonic, which needs to be taken into account 
in addition to the contribution of the zero harmonic. 

It is important to note that the scaling hypothesis is 
known to be valid at the band center. On the other hand, 
the phases 9 n are not distributed randomly in this case. 
Therefore, for energies close to E = the random phase 
approximation which is often assumed to be a core of 
the SPS conjecture, is not valid. This fact supports the 
statement of Refs. [2] that the origin of the SPS is not in 
the randomness of phases. 

A more complicated situation arises for energies close 
to the band edges, S = 2 — \E\ <C 1. However, even in this 
case the non-perturbative study (with respect to A n ) of 
classical trajectories of the map (2) allows to find the an- 
alytical expression for the Lyapunov exponent A(6) [15]. 

The HM approach is applicable also for the analysis of 
the transport properties of finite samples. The transmis- 
sion coefficient Tj, of a sample of size L can be expressed 



in terms of the radial variable as follows [19], 



l + i(r 



Li + n 2 ) 



(7) 



Here r\ i=1 2 are the radial coordinates of the last 
points [n = L) of the two trajectories that start from 

(r^,0^) - (1,0) and (r< 2) ,0< 2) ) = (1,tt/2) respec- 

tively. The values of r\ x and r L2 are calculated numer- 
ically by iterating the map (3). From equations (3) one 
gets, 



L-l 



LI 



n K< 



L-l 



n=0 



n W 



(8) 



n=0 



In what follows we apply Eq.(7) for uncorrelated ran- 
dom potential in two limit cases of ballistic and localized 
regimes. We assume that the site energies e n are dis- 
tributed randomly and homogeneously within the inter- 
val \e n \ < W/2. We consider the case of a weak disorder 
when the variance (e 2 ) = W 2 /12 is small, (e 2 ) <C 1. 
Therefore, the inverse localization length is given by, 



r 1 - 



W 2 



96(1 



E 2 \ 
4 J 



(9) 



for energies not very close to band edges and to the center 
of the energy band. 



III. BALLISTIC REGIME. 

Let us start with the ballistic regime for which the size 
L of a sample is much less than the localization length 



^ OO 



(10) 



In this case it is convenient to introduce a new parameter, 



Z L = 



Ll 



+ r 



L2 



(11) 



which is close to unity since both radii TL,i=i j 2 increase 
in "time" L very slowly, tl ~ exp(A). Therefore, the 
transmission coefficient 



l + Z L l + exp(lnZ L ) 



(12) 



can be evaluated perturbatively in terms of lnZ^. One 
can see that the statistical properties of the transmission 
coefficient Tl are entirely determined by the properties 
of Zl. The latter is related to the two trajectories of the 
classical map. In the ballistic regime a trajectory of this 
map exhibits fast rotation over angle 9 and slow diffusion 



3 



in radial direction (or, the same, in energy of the classical 
oscillator). 

In the first line we are interested in the mean value (Tl) 
of the transmission coefficient, and in its second moment 
(T|). Keeping the terms up to O (in 2 Zl) in the expan- 
sion of Tl, we obtain that the quadratic term O (in 2 Zl) 
does not contribute to Tl, 

T L «l-ilnZ L + 0(ln 3 Z L ) . (13) 

Similarily, we get, 

T 2 «l-lnZ L + iln 2 Z L + 0(ln 3 Z L ) . (14) 

After some straightforward calculations [20] that in- 
volve the map (3), the following expression for the mean 
value of In Z L is obtained, 



where 



IL-l 



(15) 



(16) 



and 



\n=0 



IL-l L-l 



\n>k k=0 



S 4 = (EE^ sin (2^)) sin (2^)) 
sin(2^)sin(2^)\ . 



(17) 



Here the terms 52 and S4 describe the correlations be- 
tween the phases 9^ and 9^ of the two classical tra- 
jectories that start from two complementary initial con- 
ditions, see Section 2. The presence of these correlation 
terms in Eq.(15) is an important fact that strongly re- 
stricts the analytical treatment. 

Let us analyze the term S2. Taking into account that 
the fluctuations of A n and 9 n are statistically indepen- 
dent, we get, 



L-l 



S2 — 81^ ^2 ^2(A„), 



(18) 



n=0 



where we introduced the two-point correlator R2 which 
depends on the scaling parameter A„ = n/lao, 



R. 



! (A„) = (sin (2^«) sin (2^))) 



(19) 



Here the average is taken over the disorder for a fixed 
number of kicks n = 1, L. 

In Fig. 1 we show numerical data for the correlator Ri 
for a wide range of the parameter A„ that covers metallic, 



A„ -c 1, and localized, A„ 3> 1 regimes. In average, the 
correlator R2 changes from —1/2 for the ballistic regime 
to 1/2 for the localized regime. This graph shows that 
the correlations give different contributions in the ballis- 
tic and localized regimes. 
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FIG. 1. Numerical data for the correlator R2 versus the 
scaling parameter A n = n/loo- The transition from the ballis- 
tic to localized regime is shown for for E = 1.5 and W = 0.1. 
The average was done over 10 4 realizations of the disorder. 
An additional "window moving" average was performed in 
order to reduce fluctuations. 
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FIG. 2. Numerical data (open circles) for fe(A n ) versus 
Eq.(21) (solid curve) for the same parameters as in Fig.l. 
and X n <C 1 



In order to evaluate analytically the correlator R2 in 
the ballistic regime, we use the approximate map for the 
angle 6 n . It is obtained from Eq.(3) in the limit |e„| <C 1, 



0k+i = 6k — H + efe 



sin 2 8k 



sin/x 



(20) 
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Using this recursion, one can express an angle 9 n in terms 
of the amplitudes eo, e...., e„_i of all the previous kicks for 
a fixed value of it. Then the following expression for i? 2 
can be obtained [20], 



i? 2 (A„) = -i+4A„-16A 2 . 



(21) 



In Fig. 2 we plot this formula together with numerical 
data for the ballistic regime, A„ <C 1. One can see that 
there is a complete agreement between numerical and 
analytical results. 

The four-point correlator S4, see Eq.(17), can be cal- 
culated analytically in a similar way [20], 



S* 4 (A) = 8A 2 + 0(A 3 ) 



(22) 



Substitution of Eqs.(21) and (22) into Eq.(15), gives the 
following formula for the mean value of In Z L 



(\nZ L ) «4A-8A 2 + 0(A 3 



(23) 



In the same way we calculate the mean value of the sec- 
ond moment, 



(ln 2 Z L ) w32A 2 + 0(A 3 ) 



(24) 



Substituting Eqs.(23) and (24) into Eqs.(13) and (14) 
respectively, we obtain, 



(T L ) = 1 - 2A + 4A 2 + O (A 3 ) 



and 



(T|) = 1 — 4A + 16 A 2 + O (A 3 ) . 
As a result, the variance reads as follows, 

Var(T L ) = (T 2 ) - (T L ) 2 = 4A 2 + O (A 3 ) 



(25) 



(26) 



(27) 



The latter expression for Var(Tj,) is known in the litera- 
ture (see, for example, [22]). However, from the analysis 
of Eqs.(25) and (26) one can obtain a more accurate ex- 
pression that depends on higher powers of A. Indeed, the 
expansions in (25) and (26) can be considered as asymp- 
totics of the following "exact" formulas, 



and 



(Tl) 



1 



1 + 2A 



1 



1 + 4A 



(28) 



(29) 



We have tested these expressions and found that they fit 
the numerical data much better than Eqs.(25) and (26) 
[20]. By combining Eq.(28) with (29), one can obtain the 
following expression for the variance of the transmission 
coefficient, 



Var(T L ) = 



4A 2 



(1 + 4A)(1 + 2A) 2 



(30) 



In Fig. 3 we compare different approximations for Var(Tt) 
with numerical data. It is clear that Eqs. (28-30) give a 
very good agreement. Note that the region of validity of 
the standard expression (27) obtained in the quadratic 
approximation is very narrow because the numerical co- 
efficients at higher terms (A 3 , A 4 , etc) in the expansion 
of Var(Tt) grow rapidly. 



0.04 




0.00 



0.00 



0.05 



0.10 



0.15 



FIG. 3. Analytical estimates of the mean variance 
Var(Ti) plotted against the numerical data (open circles) 
for E = 1.5 and W = 0.1, with the average over 
10 4 realizations of disorder. Dots stand for the standard 
quadratic approximation (27), dashed lines represent the es- 
timate that takes into account terms up to the sixth power 
of A in the expansions of (\iiZl) and (]ii 2 Zl'), namely, 
Var(T L ) = 4A 2 - 32A 3 + 176A 4 - 832A 5 + 3648A 6 . The full 
curve corresponds to our expression (30) . 



IV. LOCALIZED REGIME. 

In sufficiently long samples a regime of strong local- 
ization, A = L/loc 3> 1, is realized. Unlike the previous 
case, now all the trajectories of the classical map (3) move 
off the origin of the phase space very fast. This means 
that the radii Tl,%=\,2 are rapidly growing functions of 
discrete time L. Therefore, in this case Zl 1 and it 
is convenient to represent the transmission coefficient Tl 
in the form, 



1 



1 + Z L 



Z L 1 + ^ 



VL 



l + VL 



(31) 



Here we introduced a small parameter rjL = I/Zl <C 1 in 
order to develop a perturbative approach. In the lowest 
approximation we have, 



5 



2 " = 5T 



(32) 



It is known that in a strongly localized regime the trans- 
mission coefficient is exponentially small and reveals very 
strong fluctuations. For this reason, Tl is not a self- 
averaged quantity, and it is worth to study the logarithm 
of Tl. Its mean value is given by 



In 2- (InZ 



L 



(33) 



where the mean value of In Zl again can be expressed in 
terms of the classical trajectories of the map (3), 



(lnZ L ) =2A-ln2 + J R 



(34) 



Here the effect of correlations between the phases 8 

(2) 

and 9n enters in the last term, 
/ ( VL-l 



(i) 



-Roo = ( Ins 1 + exp 



^A l (sin2C ) -sin20«) 



.71=0 



(35) 



The detailed analysis [20] of this expression shows that, 
in fact, the term is independent of the sample length 
L. To justify this we need to take into account that in 
the limit A„ — > oo the correlator (19) shown in Fig.l ap- 
proaches 1/2, that is the mean value of sin 2 6. Therefore, 
the phases 8^ and 8^ fluctuate coherently in such a 
way that 



sin 



(2#>) - sin (20«) -> 



(36) 



for n — > oo (this result is analytically proved in [20]). 

For this reason the upper limit in the sum in Eq.(35) 
can be replaced by L = oo. Then the correlation term 
Roo becomes L— independent. The latter indicates that 
this term is Zoc-indcpcndcnt, as well, due to the scaling 
dependence of Tl on the parameter A. Being A indepen- 
dent, the correlation term R^ is a constant. Unfortu- 
nately, we are unable to evaluate this term analytically. 
The main difficulty is that the strength of the correlations 
between phases changes along the trajectory, see Fig.l. 
At the initial stage of evolution when the correlator (19) 
is different from 1/2, the two trajectories are neither sta- 
tistically independent nor coherent. This initial stage 
gives a contribution, which being small in the localized 
regime, nevertheless clearly shows that the statistics of 
\vlZl and lnTj, is not exactly Gaussian (see below). 

We evaluated the term R^ numerically and obtained 
that the following relation holds with a high accuracy, 



R^ - 2 In 2 = 



(37) 



Substituting Eq.(37) into Eq.(33) we get the standard 
expression for the mean value of lnTY, 



(lnT L ) =-2A. 



(38) 



Since this formula takes into account the correlations 
along the whole trajectory, it is accurate up to the zero 
order term with respect to A. 

In the localized regime the transmission coefficient and 
its logarithm exhibit strong fluctuations. The approx- 
imate distribution function for In Zl can be easily ob- 
tained if we neglect the initial stage of the trajectory and 



substitute r 
tain 



(i) 



„(2) ; 



in Eq.(4). Then, for In Ex we ob- 



ln Z L w t\ = ^ In D 2 n w ^ A n sin 26 n + T 2 , (39) 

n—l n—1 

where 

r 2 = J2(-^A 2 n sm 2 26 n + A 2 n sm 2 9 n SS J . (40) 

The first term in the last form of Eq.(39) is a sum of 
L >> 1 random independent numbers. Therefore, the 
statistical distribution of \uZl is the Gaussian if the 
quadratic terms T 2 are neglected. Thus, we can con- 
clude that the log-norm distribution for the transmission 
coefficient in the localized regime is obtained in the low- 
est approximation with weak disorder, and neglecting the 
difference (36) between phases along the two trajectories. 

The parameters of the Gaussian distribution for In Zl 
are calculated from Eq.(39), where the first (linear) term 
has the zero mean, however, a wide dispersion. The non- 
zero correction comes from the second (quadratic) term, 
(In Z L ) = (S 2 ) ~ \L (A 2 n ) = 2A. Calculating the second 
moment of \tiZl we can neglect the quadratic term in 
Eq.(39) and get, 



(\n 2 Z L ) «4A+(2A) 2 . 
Then, for the variance of In Zl we have, 
a 2 = (\n 2 Z L ) - (\nZ L ) 2 ^ A\, 



(41) 



(42) 



Now we can write the distribution function for In Zl , 



V (\nZ L ) = 



V27TO- 2 



exp 



(lnZ L - (lnZ L )Y 

2<7 2 



• (43) 



In Fig. 4 we fit the numerical data obtained for A = 10 
by the Gaussian distribution Eq.(43), using the disper- 
sion u as a fitting parameter. The best fit was obtained 
for a 2 = 4A + C with C = 3.2, and (\nZ L ) = 2A + In 2. 
The small correction Co to the dispersion 4A originates 
from the initial stage of the trajectories that is neglected 
in Eq.(42). At the same time, the center of the distribu- 
tion (the mean value) is in agreement with Eq.(34) where 
the calculation is performed exactly. 
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FIG. 4. Numerical data (broken curve) for the probability 
distribution of IuZl plotted against Eq.(43) (smooth curve). 
Numerical data were obtained for E = 1.5, W = 0.1 and 
A = L/loa = 10 for 10 4 realizations of the disorder. 

It should be noted that if the parameter A is not very 
big, the numerical histogram for V(\tiZl) manifests an 
asymmetry. The left tail reveals a natural cutoff at 
\uZl ~ since \uZl is a positive function, see Eqs.(7) 
and (11). This asymmetry is not visible as long as the 
dispersion a = 2\/\ is much less than the mean value 2A. 

The log-norm distribution Eq.(43) leads to the follow- 
ing mean value for the transmission coefficient [20], 




This result is slightly different from the exact formula 
(see, e.g, [21,22]), 

Having the same exponential dependence, the approxi- 
mate expression (44) differs from the exact one by an 
extra factor (tt/2) 2 (L/Zqo) - • This discrepancy origi- 
nates from the above-mentioned contribution of the ini- 
tial stage of evolution that is neglected in the distribution 
function Eq.(43). 

V. CONCLUSIONS. 

We have studied the transport properties of the ID 
standard Anderson model with weak random potential, 
using the Hamiltonian map approach. This approach is 



based on a reduction of the quantum model to the clas- 
sical two-dimensional map that describes the dynamics 
of a linear parametric oscillator with a delta-kick time 
dependence of its frequency. Amplitudes of the kicks are 
determined by the site potential of the original quantum 
model, and the energy of an eigenstate enters into the 
map as a free parameter. 

Some results have been already obtained with the use 
of this approach in application to both uncorrelated and 
correlated random potentials [7,9,10,14,19,15,16]. All 
these studies refer to the properties of the localization 
length in infinite samples. In contrast to previous re- 
sults, in this paper we consider a new question about the 
effectiveness of the Hamiltonian map approach in appli- 
cation to transport properties of finite samples. Specifi- 
cally, we are interested in the mean values of the trans- 
mission coefficient Tl, its second moment, variance, and 
the distribution functions. 

We performed an analytical treatment for two limit 
cases of ballistic and strongly localized regimes assum- 
ing weakness of the random potential. For the ballistic 
regime (when 1^ 3> L) we were able to derive analytical 
expressions for the mean values (Tl) , <T|) and the vari- 
ance (Tl) — (Tl) 2 , that are much more accurate that the 
standard estimates known in the literature. The analysis 
has revealed a non-trivial role of the correlations between 
two complimentary classical trajectories that determine 
the transmission coefficient. Our numerical study con- 
firms the analytical predictions. 

For strongly localized regime, <C L, our main inter- 
est was in the mean value of the logarithm of Tl and in 
the distribution function for lnT^. We have found that 
the leading term for (lnT^) can be easily obtained if one 
neglects the initial stage of the evolution of the comple- 
mentary classical trajectories. It is also easy to show 
that the distribution of In Tl has the Gaussian form. In 
terms of the classical map, this log-norm distribution re- 
sults from the central limit theorem applied to the ex- 
pression for the radius of classical trajectories. Corre- 
spondingly, the approach easily reproduces the estimates 
for the mean value of the second moment of InTj, and 
for its variance. 

The results of our analysis may find further applica- 
tions to the problem of the single parameter scaling, as 
well as in the study of transport properties of ID random 
models with long-range correlations in the site potential. 
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